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ABSTRACT 

We present a new and simple method to measure the instantaneous mass and radial 
growth rates of the stellar discs of spiral galaxies, based on their star formation rate 
surface density (SFRD) profiles. Under the hypothesis that discs are exponential with 
time-varying scalelengths, we derive a universal theoretical profile for the SFRD, with 
a linear dependence on two parameters: the specific mass growth rate um = M*/M* 
and the specific radial growth rate zzr = of the disc. We test our theory on a 

sample of 35 nearby spiral galaxies, for which we derive a measurement of vm and zzr. 
32/35 galaxies show the signature of ongoing inside-out growth (zzr > 0). The typical 
derived e-folding timescales for mass and radial growth in our sample are ~ 10 Gyr 
and ~ 30 Gyr, respectively, with some systematic uncertainties. More massive discs 
have a larger scatter in I'm and z/r, biased towards a slower growth, both in mass 
and size. We find a linear relation between the two growth rates, indicating that our 
galaxy discs grow in size at ~ 0.35 times the rate at which they grow in mass; this 
ratio is largely unaffected by systematics. Our results are in very good agreement with 
theoretical expectations if known scaling relations of disc galaxies are not evolving 
with time. 
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1 INTRODUCTION 


(Freeman 1970), though sometimes broken at galaxy edges 


The theory of cosmological tidal torques (Peebles|l969) pre¬ 
dicts the mean specific angular momentum of galaxies to be 
an increasing function of time. If applied to spiral galaxies, in 
which stars are mostly distributed on a rotating, centrifu- 
gally supported, disc, the theory suggests that the outer 
parts, with higher specific angular momenta, should form 
later than the inner ones ( Larson|1976| the so-called inside- 
out formation scenario), also implying that spirals should 
grow in size while they grow in mass. Apart from this quite 
general prediction provided by cosmology, the details about 
how stellar discs form and grow in mass and size are not 
known from first principles and significant observational ef¬ 
fort is still required to shed light on the missing links from 
structure formation to galaxy formation. 

An invaluable input for modelers comes from some sim¬ 
ple observed properties of the discs of spiral galaxies, that 
still wait for a comparatively simple theoretical explanation: 
among them, the exponential radial structure of galaxy discs 


* E-mail: gabriele.pezzulli@unibo.it 


(e.g. Pohlen & Trujillo 2006 Erwin, Pohlen & Beckman 


2008), and the fact that they obey simple scaling relations, 
including the Tully-Fisher relation between rotational ve¬ 
locity and mass ( Tully fc Fisher|1977 see also the ‘baryonic 
Tully-Fisher relation’, |McGaugh|2012||Zaritsky et al.|2014[), 


the Fall relation between angular momentum and mass (Fall 
1983 Romanowsky fe Fall|2012 l and a more scattered mass- 


size relation (e.g. |Shen et al.||2003| |Courteau et al.||2007| ), 
which can also be seen as a corollary of the other two. 

Observations of galaxies at different redshifts indicate 
that stellar discs have an exponential structure since very 
early epochs (Elmegreen et al. 2005; Fa thi et al.|2012 1, while 
it is less clear whether scaling relations are truly universal or 
they evolve with cosmic time. For example, direct measure¬ 
ments of the mass-size relation for disc galaxies at various 
redshifts has led to claims for little or no evolution (e.g. 


Ravindranath et al. 2004; Barden et al. 2005 ; Ichikawa, Ka- 

jisawa & Akhlaghi 2012 

1 as well as significant or strong evo- 

lution (e.g. 

Mao, Mo & White 1998 Buitrago et al. 

2008 

Fathi et al. 

2012). The interpretation and comparison of 
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2 Pezzulli et al. 


geneities among observations at different redshifts, as well as 
differences in sample definitions and analysis techniques (see 
e.g. Lange et al. 20151; also, several subtle issues have been 
shown to significantly bias the results, most notably the se¬ 


lection effect due to cosmological dimming (e.g. Simard et al. 


1999) and the evolution of M/L ratios due to evolving stellar 
populations (e.g. Trujillo et al.|[2006). The related problem 


of the possible evolution of the Tully-Fisher relation, which 
also involves kinematic measurements, is even more complex 


and controversial (e.g Vogt et al.||1997| Mao, Mo & White 
|1998| |Miller et al.|[2011| |Miller et al.||2012| and references 
therein). 

Since these issues are of extreme importance for our 
understanding of galaxy evolution, more observational ef¬ 
fort is desirable, possibly dealing with multiple independent 
probes, to unveil the growth of the exponential discs of spi¬ 
ral galaxies. In addition to the crucial, but often challeng¬ 
ing, comparison of galaxy properties at different redshifts, 
indirect information can be gained on the size growth of 
galaxy discs from the study of their properties in the Local 
Universe. Efforts in this direction can be split in two cat¬ 
egories. The first is the exploitation of fossil signals of the 
past structure of the disc: most notably, chemical enrichment 


(e.g. Boissier & Prantzos| 

1999 Chiappini, Matteucci & Ro- 

mano][2001 Molla & Diaz 

2005 Naab & Ostriker 2006), and 


properties of stellar populations, including colour gradients 


(e.g. Bell & de Jong 2000 MacArthur et al. 2004 Wang 
et al. |2011 1, spectrophotometry (e.g. Munoz-Mateos et al. 


2011 |Gonzalez Delgado et al. 


diagrams (e.g. Williams et al.||2009 Gogarten et al. 2010 


and colour-magnitude 


Barker et al. 20111. The second possibility is to look for 


the instantaneous signal of the growth process while it is 
in act. Spiral galaxies are not just passively evolving stellar 
systems, but keep forming stars at a sustained rate through¬ 
out their evolution (e.g. Aumer & Binney 2009 Fraternali & 


Tomassetti| f2C) 12 Toje iro et al.| |20131. Therefore, the radial 

distribution of newly born stars is a crucial ingredient for 
the structural evolution of a stellar disc and it can be used 
as a clean and direct probe of its growth. 

Thanks to the deep UV photometry of the Galaxy Evo¬ 
lution Explorer (GALEX, M artin et al.|2 005 radial profiles 
of the star formation rate surface density (SFRD) of nearby 
galaxies can now be traced out to considerable galactocen- 
tric distances and low levels of star formation activity. SFRD 
profiles, as traced by the UV light emitted by young stars, of¬ 
ten turn out to be quite regular and, in many cases, remark¬ 


ably similar to exponentials (Boissier et al. 2007 see also 


Goddard, Kennicutt & Ryan-Weber 2010). This supports 


the idea that star formation is tightly linked with whatever 
process is responsible for the mantainance and evolution of 
the exponential structure of galaxy discs. A closer inspec¬ 
tion of the aforementioned SFRD profiles reveals, in many 
galaxies, some deviations from the exponential shape, in the 


form of a central downbending or depletion (see also Munoz- 
Mateos et al. 2011). This is also clearly seen in the SFRD of 
the Milky Way, as traced, for example, by the distribution of 
Galactic supernova remnants (Case & Bhattacharya 19981. 
Indeed, star formation becoming progressively less effective 
towards the centre of galaxies is not very surprising, within 
the inside-out formation scenario: in the inner regions, the 
bulk of gas accretion and conversion into stars occurs quite 
early and relatively little residual star formation is expected 


to be in place there at late epochs, while the outskirts are 
still in their youth. Ultimately, the observed properties of 
SFRD profiles of spiral galaxies are in qualitative agreement 
with the inside-out paradigm and therefore they are good 
candidates to enclose the signal of radial growth. The aim 
of this work is to give a simple quantitative description of 
this signal and a method to derive a measurement of the 
instantaneous mass and radial growth rates of the discs of 
spiral galaxies from the amplitude and shape of their SFRD 

this direction has been done by 
. They assumed that surface den- 
and star formation rate can be 
approximated with exponential profiles, though with differ¬ 
ent scalelengths, implying exponential profiles for the spe¬ 
cific star formation rate (sSFR) as well, with sSFR in¬ 
creasing with radius for inside-out growing galaxies. This 
parametrization was applied to a sample of nearby spiral 
galaxies and the results were compared with the predictions 
of simple structural evolution models, providing constraints 
on the inside-out process. 

In this work we make a step forward, proposing a 
method that is both simpler and more powerful. Rather 
than modeling both stellar mass and SFRD with exponen¬ 
tials, we assume that just stellar discs are exponential at any 
time, with time-varying scalelengths. This naturally brings 
us (Sec. [2| to predict a universal shape for the SFRD pro¬ 
file with the observed properties outlined above, namely an 
inner depletion and an outer exponential decline. Further¬ 
more, our theoretical profile has a very simple (linear) de¬ 
pendence on the disc mass and radial growth rates; hence, 
these parameters can be directly derived from observations 
in a model-independent way. We apply our method to a sam¬ 
ple of nearby spiral galaxies described in Sec. [ 3 ] discuss our 
analysis in Sec. [4] and present our results in Sec. [5] In Sec. 
[6] the consequences of our findings are discussed on the is¬ 
sue of whether known scaling relations for galaxy discs are 
evolving with time or not. In Sec. [T] we give a summary. 


profiles. 

An earlier attempt in 


Munoz-Mateos et al. (2007 


sities of both stellar mass 


2 STAR FORMATION IN EXPONENTIAL 
DISCS 

2.1 A simple model 


Let us assume that the mass surface density E* of the stellar 
disc of a spiral galaxy is well described, at any time, by an 
exponentially declining function of radius R, identified by a 
radial scalelength R* and a mass M* Q both allowed to vary 
with time t: 


E* {t , R) 


M*(t) „( R \ 

27tR*(£) 6XP \ R*(t)J 


(1) 


Just taking the partial time derivative of |l]) we get a very 
simple prediction for the star formation rate surface density 
E* as a function of time and galactocentric radius: 


E* (t, R) 


(^m (t) + m(t) 



S* (t, i?) 


( 2 ) 


1 Throughout the paper, when referring to a spiral galaxy, we 
will use the symbol M* to denote the stellar mass of its disc 
component alone. 
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Radial growth of stellar discs 3 


where E* is given by 0 , while the quantities um and ur are 
defined by: 


/ \ d /, , , , .. MJt ) 
u m {t) :=_(lnM*(t)) = ^ 


Mt) :=^(JnA(*)) = ^|| 


We discuss them more thoroughly in Sec. |2.3| 


(3) 

(4) 


effects can sometimes be observationally inferred in low sur¬ 
face brightness regions or in low surface brightness galaxies 
as a whole (Boissier et al.|2008l and may also be related to 


the phenomenology of extended UV (XUV) discs (Thilker 
et al. [2007a I. Unfortunately, in general, neither the magni¬ 


tude nor the direction of the needed correction can be known 
a priori. However, these possibilities should be kept in mind 
when considering peculiar features in the observed SFRD 
profiles of some individual objects. 


2.2 Theoretical caveats 

At least two caveats should be kept in mind when consider¬ 
ing the elementary inference outlined in Sec. |2.1| 

First, by identifying E* with (<?E */dt), we have implic¬ 
itly neglected any contribution coming from a possible net 
radial flux of stars, that is l/R d/8R{2-KRT, i ,u I f), u+ being 
the net radial velocity of stars. While radial migration of 
stars is widely recognized to be a fundamental ingredient of 


galaxy evolution, it has also been shown (Sellwood & Binney 
2002 Roskar et al.|2012 ) that its main working mechanism is 


basically a switch in the radial position of two stars in differ¬ 


ent circular orbits (the so-called churning , Sehonrich & Bin- 
ney||2009 1. This process produces no dynamical heating, no 


net radial flow of stars and no change in the mass distribu¬ 
tion of the disc. Of course, some minor contribution to radial 
migration are also expected from other processes: breaks to 
our approximation can be expected in some cases, mostly in 
the inner regions, where dynamical processes might be as¬ 
sociated with the formation of bars, rings and pseudobulges 
(e.g. Sellwood| 20141 and at the outer edge, where radial mi¬ 
gration has been proposed to induce changes in the outer 


structure of discs (Yoachim, Roskar, & Debattista 2012). 


More complex effects are also possible due to the interplay 
between stellar dynamics, gas dynamics and star formation; 
quite different approaches to this problem can be found, for 


example, in Sehonrich & Binney (2009), Kubryk, Prantzos 


(2014). 


& Athanassoula (20131 and Minchev, Chiappini & Martig 


Second, since stellar populations, during their evolution, 


return a substantial fraction of their mass to the ISM (Tins- 
ley][l980), it is not necessarily trivial to connect the time 
derivative E* to observed values of SFRD. In the follow¬ 
ing, we will adopt the instantaneous recycling approximation 
(IR A) and assume that a constant return fraction 1Z of the 
mass of a stellar population is instantaneously given back 
to the ISM. Accordingly, our E* represents the reduced star 
formation rate surface density and it is equal to the observed 
SFRD multiplied by a factor (1 — 1Z), although we will often 
omit the attribute reduced , for brevity. More detailed studies 
(e.g. |Leitner fe Kravtsov|2011| ) show that the majority of the 
returned mass is released within ~ 1 Gyr from the birth of a 
population; hence, our approximation will be valid in those 
galaxies, or galaxy regions, where star formation has not 
been varying abruptly on timescales shorter than ~ 1 Gyr. 
It can be easily seen that such abrupt changes can, in prin¬ 
ciple, be taken into account by replacing 7 Z with an effec¬ 
tive return fraction 7Z e s, which is higher or lower than TZ for 
abrupt quenching or starbursting events, respectively. These 
may be due, for instance, to tidal or ram pressure stripping 
or, viceversa, to significant recent accretion events. Similar 


2.3 Mass and radial growth rates 

The quantities um and ur, defined in <0 and Q, shall be 
called specific mass growth rate and specific radial growth 
rate. The word specific refers to the fact that they represent 
the mass and radial growth rates M* and 7?.*, normalized to 
the actual value of mass M* and scalelength f?*, respectively. 
However, since in this work we will deal only with specific 
quantities, we will often omit the attribute and refer to them 
just as mass and radial growth rates, for brevity. 

While um is always positive (stellar mass is never de¬ 
stroyed), vr can in principle take both signs, positive values 
being expected in the case of inside-out growth. At any time, 
the inverse of um and ur can be interpreted as instantaneous 
estimates of the timescales for the growth of the stellar mass 
and scalelength, respectively (or for disc shrinking, in the 
case ur < 0). 

We notice that um is strictly related to another physi¬ 
cal quantity, namely the specific star formation rate (sSFR). 
More precisely, following the terminology of |Lilly et al.| 
(20131, um coincides with the reduced specific star forma¬ 


tion rate of the disc, where the word reduced (which we will 
omit from now on) refers to the fact that we are including 
the effect of the return fraction TZ. Since our um refers to 
the disc alone, it should not be confused with the sSFR of 
a whole spiral galaxy, which is evaluated including also the 
other stellar components, like the bulge. While the bulge can 
give a non-negligible contribution to the total stellar mass, 
it usually harbours little or no star formation: hence, the 
sSFR of a whole galaxy and of its disc alone can differ sig¬ 
nificantly (Abramson et al. |2014 1. Also, um should not be 
confused with the local sSFR (E*/E*), which is, in general, 


a function of galactocentric radius (e.g. Munoz-Mateos et al. 
2007). 


The analogous quantity for the stellar scalelength, the 
(specific) radial growth rate Ur, has been studied much less 
(and, as far as we know, not even clearly defined until now). 
To provide a simple method for its measurement is the main 
aim of this work. Since um and Ur have the same physical 
dimensions and refer to the two basic properties of an expo¬ 
nential disc, we are also interested to measure both quanti¬ 
ties at the same time and to attempt a comparison between 
them. This is indeed a natural outcome of our method (Sec. 
2.4 1 and will bring us to the most interesting consequences 
of our results (Sec. [6j. 


2.4 Predicted properties of SFRD profiles 

Our simple model predicts that, if a galaxy is observed at 
some particular time, its SFRD should follow a radial profile 


© 2015 RAS, MNRAS 000, [T] ?? 



















































4 Pezzulli et al. 



R/R* 


Figure 1. Predicted shape of the SFRD profile, in dimension¬ 
less units, for some illustrative values of the radial growth rate 
z/r: absence of radial growth (z/r = 0, dotted black line), inside- 
out growth (z/r = 1/3 I'm, solid red line) and disc shrinking 
(z/r = —1/3 I'm, dashed blue line). The radius is in units of i?*, 
the SFRD is normalized to Ei = z/mM*/ 27 tR*, so that the com¬ 
parison refers to discs with the same stellar mass, scalelength and 
global sSFR. 


of the form: 

+ - 2 )) exp (-f) (5) 

In Fig.[T]the predicted shape of the SFRD profile is drawn 
out, in dimensionless units, for some representative situa¬ 
tions, which differ for the sign of the radial growth param¬ 
eter j/r. To achieve a fully dimensionless description of the 
model, we use here, as a parameter, the dimensionless ratio 
i/r/j/m, which has the same sign of i/R, since, as pointed out 
in Sec. |2.3| vm is always positive. Also, with our adopted 
dimensionless units, the curves in Fig. [T] compare with each 
other as model discs that share the same mass, scalelength 
and global sSFR vm, but differ in the spatial distribution of 
star formation, depending on the presence (and the direc¬ 
tion) of an evolution of the scalelength with time. 

In the absence of radial evolution (;zr = 0), the scale- 
length of the stellar disc is constant with time, stars always 
form with the same spatial distribution and the SFRD pro¬ 
file is an exponential as well, with the same scalelength of 
the already present stellar disc. In the case of disc shrink¬ 
ing (i/R < 0) star formation is enhanced, with respect to 
the previous case, in the inner regions, but it abruptly 
drops in the outskirts (reaching a vertical asymptote at 
R/R+ = (2 + ^m/|^r|)). Conversely, for inside-out growth 
(i'r > 0), the SFRD shows the characteristic depletion in 
the central regions, while it is enhanced in the outskirts, 
where it ultimately gently declines with increasing radius, 
with an asymptotic behaviour, at large radii, similar to the 
one of the stellar mass distribution. 


in the context of the SINGS survey (Kennicutt et al.(2003|. 
For these and other nearby galaxies, radial profiles have been 
derived and published by Munoz-Mateos et al. (2009b I and 


Munoz-Mateos et al. (2009a), for multi-wavelength emission 


ranging, for most galaxies, from FIR to FUV. Such a broad 
range is useful to trace both the stellar mass and the star 
formation rate, corrected for the effect of dust extinction. 


Munoz-Mateos et al. (2009a I give radial profiles for the ex¬ 


tinction in the UV, as inferred from TIR/UV flux ratios. 

For this work, we made use of the radial profiles of the 
emission in the FUV GALEX band, corrected for extinction 
in the FUV (Apuv)- We also used profiles of emission in the 
3.6 /im IRAC band, which we assume to be a good tracer 
of the stellar mass surface density. Some contamination may 
arise from the 3.3 /rm PAH line, hot dust and AGB stars; 
however, these contributions are only expected to be impor¬ 
tant at the small scales of individual star-forming regions 


and just a mild effect persists at larger scales (Meidt et al. 


2012). Furthermore, dust extinction at this wavelength is 


negligible and the M/L is quite insensitive to variations in 


age and metallicity, if compared to the optical bands (Meidt 
et al.|20l4|. 


From the original sample of 42 galaxies, we excluded 
5 galaxies (NGC 3049, NGC 3938, NGC 4254, NGC 4321, 
NGC 4450) for which FUV meaurements were lacking be¬ 
cause the GALEX FUV detector was turned off for technical 
reasons, one galaxy (NGC 7552) which was not present in 
the Munoz-Mateos et al. (2009a) sample and another one 
(NGC 4625) for which the Afuv radial profile was consti¬ 
tuted of only one point. Therefore we ended up with a final 
sample of 35 galaxies. We have considered some possible 
additional criteria to further restrict our sample. Since our 
main goal is the study of the slow, continuous, evolution of 
the regular exponential structure of discs, galaxies that are 
suspected to be undergoing violent transient events, like in¬ 
teractions or mergers, could be excluded from the analysis. 
Furthermore, since we make a quantitative analysis of az- 
imuthally averaged radial profiles, we could exclude those 
galaxies for which the geometrical parameters involved in 
the average (inclination and position angle) are not known 
with good accuracy or are suspected to be varying with ra¬ 
dius. At least 3 galaxies (NGC 1097, NGC 1512, NGC 5194) 
have a nearby companion and for at least one (NGC 1512) 
the adopted position angle reproduces the outer isophotes 
better than the inner ones. However, it is not clear whether 
such selections could be done in a completely unbiased way. 
Also considering that our sample is relatively small, we de¬ 
cided to homogeneously analyse the whole set of 35 galaxies. 
Nonetheless, the aforementioned caveats should be kept in 
mind while considering our results. We adopt morphological 
classifications, distances and inclinations (as derived from 
axis ratios) from Munoz-Mateos et al. (2011); these proper¬ 
ties can also be found here in Tabled] 


3 SAMPLE AND DATA 
3.1 Sample definition 

To define our sample, we started from the one studied by 
Munoz-Mateos et al. (2011), which consists of 42 nearby 


spiral galaxies observed both by Spitzer and by GALEX, 


3.2 Stellar mass surface density 

To get the stellar mass distribution, we made use of the high 


resolution (6 arcsec) radial profiles at 3.6 p m from Munoz- 
|Mateos et al.| ( |2009b| ). To convert from surface brightness to 
mass surface density, we used the conversion formula: 

X* , „ „7 . I3.6iJ.u1 


M© pc - 


= 1.9 • 10' cos* 


Jy arcsec - 


( 6 ) 
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Table 1. Basic properties of galaxies in our sample. Morphological classification, distances and inclinations (derived from axis ratios) 
are as in Munoz-Mateos et al. (2011). 


Galaxy 

RA (J2000) 
h m s 

Dec (J2000) 

O / // 

Morphological 

Type 

Type 

T 

Distance 

(Mpc) 

cos i 

NGC 0024 

00 09 56.5 

-24 57 47.3 

SA(s)c 

5 

8.2 

0.224 

NGC 0337 

00 59 50.1 

-07 34 40.7 

SB(s)d 

7 

25 

0.621 

NGC 0628 

01 36 41.8 

15 47 00.5 

SA(s)c 

5 

11 

0.905 

NGC 0925 

02 27 16.9 

33 34 45.0 

SAB(s)d 

7 

9.3 

0.562 

NGC 1097 

02 46 19.1 

-30 16 29.7 

SB(s)b 

3 

15 

0.677 

NGC 1512 

04 03 54.3 

-43 20 55.9 

SB(r)a 

1 

10 

0.629 

NGC 1566 

04 20 00.4 

-54 56 16.1 

SAB(s)bc 

4 

17 

0.795 

NGC 2403 

07 36 51.4 

65 36 09.2 

SAB(s)cd 

6 

3.2 

0.562 

NGC 2841 

09 22 02.6 

50 58 35.5 

SA(r)b 

3 

14 

0.432 

NGC 2976 

09 47 15.5 

67 54 59.0 

SAc pec 

5 

3.6 

0.458 

NGC 3031 

09 55 33.2 

69 03 55.1 

SA(s)ab 

2 

3.6 

0.524 

NGC 3184 

10 18 17.0 

41 25 28.0 

SAB(rs)cd 

6 

8.6 

0.932 

NGC 3198 

10 19 54.9 

45 32 59.0 

SB(rs)c 

5 

17 

0.388 

IC 2574 

10 28 23.5 

68 24 43.7 

SAB(s)m 

9 

4.0 

0.409 

NGC 3351 

10 43 57.7 

11 42 13.0 

SB(r)b 

3 

12 

0.676 

NGC 3521 

11 05 48.6 

-00 02 09.1 

SAB(rs)bc 

4 

9.0 

0.464 

NGC 3621 

11 18 16.5 

-32 48 50.6 

SA(s)d 

7 

8.3 

0.577 

NGC 3627 

11 20 15.0 

12 59 29.6 

SAB(s)b 

3 

9.1 

0.462 

NGC 4236 

12 16 42.1 

69 27 45.3 

SB(s)dm 

8 

4.5 

0.329 

NGC 4536 

12 34 27.1 

02 11 16.4 

SAB(rs)bc 

4 

15 

0.421 

NGC 4559 

12 35 57.7 

27 57 35.1 

SAB(rs)cd 

6 

17 

0.411 

NGC 4569 

12 36 49.8 

13 09 46.3 

SAB(rs)ab 

2 

17 

0.463 

NGC 4579 

12 37 43.6 

11 49 05.1 

SAB(rs)b 

3 

17 

0.797 

NGC 4725 

12 50 26.6 

25 30 02.7 

SAB(r)ab pec 

2 

17 

0.710 

NGC 4736 

12 50 53.1 

41 07 13.6 

(R)SA(r)ab 

2 

5.2 

0.813 

NGC 4826 

12 56 43.8 

21 40 51.9 

(R)SA(rs)ab 

2 

7.5 

0.540 

NGC 5033 

13 13 27.5 

36 35 38.0 

SA(s)c 

5 

13 

0.467 

NGC 5055 

13 15 49.3 

42 01 45.4 

SA(rs)bc 

4 

8.2 

0.571 

NGC 5194 

13 29 52.7 

47 11 42.6 

SA(s)bc pec 

4 

8.4 

0.804 

NGC 5398“ 

14 01 21.6 

-33 03 49.6 

(R’)SB(s)dm pec 

8.1 

16 

0.607 

NGC 5713 

14 40 11.5 

-00 17 21.2 

SAB(rs)bc pec 

4 

27 

0.893 

IC 4710 

18 28 38.0 

-66 58 56.0 

SB(s)m 

9 

8.5 

0.778 

NGC 6946 

20 34 52.3 

60 09 14.2 

SAB(rs)cd 

6 

5.5 

0.852 

NGC 7331 

22 37 04.1 

34 24 56.3 

SA(s)b 

3 

15 

0.352 

NGC 7793 

23 57 49.8 

-32 35 27.7 

SA(s)d 

7 

3.9 

0.677 


(a) In the original sample, this galaxy was referred to as TOL 89, which is the name of an HII region embedded within it. 



3.3 Star formation rate surface density 

To derive the SFRD profiles, we took the low resolution (48 


arcsec) radial profiles in the FUV band from Munoz-Mateos 
et al. (|2009a I and corrected them for extinction using the 


Afuv radial profiles at the same resolution. In that work, 
two estimates of Afuv are provided, based on two slightly 


2 The formula of |Leroy et al.| j2008| was based on their measured 
linear 3.6 pm-to-K-band flux conversion and on an assumed K- 
band mass-to-light ratio. Changing the latter from 0.5 to 0.8 is 
equivalent to introducing an additional factor 1.6 in the 3.6 /im- 
to-stellar mass conversion. 


different dust attenuation prescriptions by Buat et al. (20051 


and by Cortese et al. (20081, the latter containing a refine¬ 


ment to take additional dust heating from old stars into ac¬ 
count. In this work we used extinction profiles from the [Buat] 
et al. (2005) prescription. We made this choice to maximize 


simplicity and reproducibility of our analysis (this recipe 
does not require additional information on K-band photom¬ 


etry). Munoz-Mateos et al. (2009a) showed that the two pre¬ 


scriptions differ significantly only for early-type galaxies (el¬ 
lipticals and lenticulars), which are absent in our sample.We 
verified that our conclusions are not modified when changing 
the adopted prescription. 

The extinction corrected profile /xpuv.corr (expressed in 
the AB magnitude system) was converted into a SFRD by 
making use of the formula: 


— = (1 — 72.) cos i 10 


-0.4/ipuV,corr + 10.413 


(7) 


M© pc -2 Gyr 

which is again consistent with Munoz-Mateos et al. (20071. 


For the return fraction, we adopted 72. = 0.3, which is inter¬ 
mediate between possible values for different IMF choices 


(Leitner & Kravtsov 2011 Fraternali & Tomassetti 2012). 


© 2015 RAS, MNRAS 000, [T] ?? 












































6 Pezzulli et al. 


NGC 0628 
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Figure 2. Our analysis for the galaxy NGC 628. Upper panels Selection of the domain where the emission at 3.6 pm is dominated by the 
light from the stars in the the disc: left the inner ellipse at 33 arcsec, out of which the spiral structure appears, right the outer ellipse at 
261 arcsec, out of which the contribution from noise becomes significant. Note that the two images have very different scale and contrast. 
The inner ellipse is also shown in the right panel, to make the whole selected region visible at once. Lower-left panel Exponential fit to 
the stellar mass surface density, as traced by the emission at 3.6 ^m. The vertical dashed lines mark the limits of the domain that we 
have selected for this fit; in this case, the best-fitting exponential also extends further out in the outskirts. Lower-right panel Fit of our 
theoretical SFRD profile to the one obtained from extinction-corrected FUV light; a visual comparison with Fig. ^ is already enough to 
recognize this as an inside-out growing galaxy. 


Also the systematics associated with 0 is discussed in Sec. 

EGO 


ters: the specific mass growth rate um and the specific radial 
growth rate !/r. 


4 ANALYSIS 

For each galaxy in the sample, we performed our analysis 
in two steps. First, we made an exponential fit to the radial 
profile of the stellar mass surface density of the disc (see 
Sec. 3.21, deriving the values for the disc mass M* and scale- 
length R+. Then, keeping these parameters fixed, we fitted 
our theoretical profile 0 to the SFRD data (see Sec. |3.3| ). 
This second fit is the test bed for our theory. If successful, it 
provides our measurement of the two disc growth parame- 


A more extended description of the two steps is given 
in Secc. |4.1| and |4.2| They are depicted, for each individual 
galaxy, in an Atlas, which we provide as supplementary on¬ 
line material. A representative example, for the galaxy NGC 
628, is reported here for illustrative purposes (Fig. 0. 


All our fits were performed with a standard Marquardt- 
Levenberg algorithm. We also repeated the whole analysis 
with a different method (see Sec. 4.31 and verified that our 
results are robust with respect to the fitting strategy. 
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4.1 Fit of exponential discs 


In order to extract the disc parameters from the radial pro¬ 
files described in Sec. |3.2[ we performed a simple exponential 
fit, for each galaxy, on a radial domain where the NIR emis¬ 
sion is dominated by the disc component. 


Such a domain was identified, on a case-by-case basis, 
considering the shape of the 3.6 /im profile with the aid of 
the direct visual inspection of the 2D maps at the same wave¬ 
length. For details about how these maps were obtained the 
rea der is referred tojRegan et al~| ( |2004| ), |Dale et al.| ( |2005| ) 
and Munoz-Mateos et al. (2009b I. For each galaxy, the min¬ 
imum and the maximum radius of our selected domain are 
the semi-major axes of two concentric ellipses, with cen¬ 
tre and orientation equal to the ones used in the derivation 
of the profiles. The inner ellipse was chosen to exclude the 
central bright component, if present, like a bulge, a bar, or 
a central ring; the detection of spiral arms has been used 
in some cases as an evidence for the prominence of the disc 
component in a given region. In 4 cases (NGC 2403, IC 2574, 
NGC 4236, NGC 4826) we have found that the adopted 
centre of the ellipses did not coincide with the peak of the 
3.6 /im emission; for these galaxies, an inner ellipse was se¬ 
lected with a semi-major axis equal or greater than twice 
the observed offset. The outer ellipse was most of the times 
selected to exclude those external regions where a contribu¬ 
tion to the emission coming from the noise was found to be 
significant; for our data, this happens at a typical value of 
log(/ 3.6 Iim/Jy arcsec -2 )-6.5. 

In 5 galaxies (NGC 3521, NGC 3621, NGC 4736, NGC 
5055, NGC 7331) we found a significant flattening of the 3.6 
/rm profile well above the noise level and we excluded the 
outer region from the exponential fit for these objects. The 
most striking case is NGC 4736, where the effect is probably 
related to the presence of a prominent outer ring. In 3 of the 
above cases (NGC 3521, NGC 3621, NGC 5055) the change 
of slope occurs very near to the outermost radius where a 
regular spiral pattern can be seen. In the remaining two 
objects (NGC 3521 and NGC 7331) the flattening is associ¬ 
ated with an abrupt change in the geometry of the 3.6 /jm 
emission, with the isophotes becoming remarkably large and 
irregular in the outer regions. We ignore the physical origin 
of this effect; nonetheless, these peculiarities should be kept 
in mind in the interpretation of our results for these objects. 


In most cases, our fits were performed weighting each 
point according to the nominal error, quoted in the origi¬ 
nal profiles. However, this is not necessarily always the best 
choice. Real galaxies are not expected to precisely follow an 
exponential, since transient perturbations like spiral arms, 
which are by definition ubiquitous in spiral galaxies, can 
sometimes overimpose oscillations on an underlying regu¬ 
lar disc. This effect can become particularly important in 
the presence of spiral arms with a small pitch angle, tend¬ 
ing to dominate the emission in a limited radial range. If, 
for a given galaxy, points with small error bars happen to 
be preferentially located in a region dominated by spiral 
structure, the formal best-fitting profile will be biased to 
reproduce transient features, potentially missing the over¬ 
all structure of the disc. For 7 galaxies (NGC 1097, NGC 
1512, NGC 3031, NGC 3184, NGC 3351, NGC 4569, NGC 
7793) we have found that an unweighted fit provided a bet¬ 


ter description of the overall structure of the profiles in the 
considered radial range. 

Although the whole procedure is slightly subjective, we 
verified that it gave a better account to the observed prop¬ 
erties of our galaxies, in the domain of interest, with respect 
to a more complex global analysis, involving more compo¬ 
nents and parameters. This approach is the most suitable for 
our purposes, since we are just interested to reliably derive 
the disc parameters, rather than to get a detailed struc¬ 
tural decomposition of the whole galaxy. We verified that 
the parameters we found were stable with respect to small 
variations in the selection of the inner and outer radii. 

An example of our domain selection and disc-fitting pro¬ 
cedure is given in Fig. [2] (upper panels and lower-left panel) 
for the case of NGC 628; similar images and plots for the 
other galaxies can be found in the online Atlas. 


4.2 Fit of the star formation rate surface density 

In the second step of our analysis, we fitted equation § 
to the observed SFRD profiles, keeping fixed the structural 
parameters and R t found in the previous step. 

The SFRD profiles have a worst spatial resolution, and 
hence a more limited number of independent points, with re¬ 
spect to the mass surface density profiles. As a consequence, 
the results of the SFRD fits are more sensitive to changes 
in the adopted radial domain. In order to limit the depen¬ 
dence of our analysis on subjective choices, we decided to 
always perform the fit on the whole available domain. Not 
to exclude any point from the inner regions is equivalent 
to assume that the bulk of star formation is everywhere as¬ 
sociated with the disc component. In other words, we ne¬ 
glected possible star formation activity directly occurring in 
the bulges, which is quite reasonable since these structures 
are known to be dominated by old stellar populations. Nei¬ 
ther we put outer limits to our domain, implying that we did 
not try to model possible transient star formation episodes 
that might dominate the UV emission in the outer regions, 
nor any kind of structural irregularity and, most noticeably, 
the possible presence of warps, which are generally expected 
in the periphery of discs (Briggs 1990). All the SFRD fits 
were performed weighting points with their nominal errors, 
which were derived just propagating the errors in the /xfuv 
and Afuv profiles. 

In Fig. [2] (lower-right panel) the best-fitting SFRD pro¬ 
file is reported for the galaxy NGC 628; similar plots are 
reported for all galaxies in the online Atlas. 

In considering this part of the analysis, it should be 
kept in mind that, while the parameters um and i/r, in the 
theoretical SFRD profile |5|, are allowed to change in the fit¬ 
ting process, the global slope is strongly constrained by the 
parameter i?*, which is held fixed to the value previously 
obtained from the structural fit (Sec. 4.1 1 . Hence, notwith¬ 
standing the presence of 2 free parameters, we are by no 
means able to reproduce arbitrary profiles and the fact that 
we can recover the majority of SFRD distributions shall be 
regarded as a success of the model and gives us confidence 
on the meaningfulness of the resulting best-fit parameters. 

In 5 cases (NGC 1512, NGC 3521, NGC 3621, NGC 
4736, NGC 5055) we clearly detect an outer flattening in 
the radial profiles of both SFRD and stellar mass surface 
density. This can be considered as an indication for the ex- 
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istence of a distinct, relatively long-lived, outer star form¬ 
ing component. More detailed studies would be necessary 
to clarify this point. However, we notice here that three of 
these galaxies (NGC 1512, NGC 3621, NGC 5055) have been 


classified by Thilker et al. (2007b) as having a Type 1 XUV 


disc. We also notice that sometimes (e.g. for NGC 3621) 
the spatial coincidence between the two breaks is perfect, 
while in other cases (NGC 3521 and NGC 5055) the SFRD 
flattening occurs at larger radii, maybe challenging the idea 
of a common origin of the two phenomena. A unique case is 
the one of NGC 7331, which has a very prominent flattening 
of the stellar mass distribution, but an almost exponential 
SFRD profile, which our model is unable to account for; we 


refer the reader to Thilker et al. (2007b) and Ludwig et al. 


(20121 for more specific studies on this peculiar object and 


its surroundings. Finally, we report one case (IC 4710) where 
a quite marked downbending is found, at the same radius, 
in both profiles. This is not a very secure result, since the 


break occurs out of our chosen outer ellipse (see Sec. 4.11; if 
confirmed, it may be an example of radial migration in the 
presence of an outer cut-off in star formation efficiency, as 


also Sec. 2.2 1 . 


described e.g. by Yoachim, Roskar & Debattista (2012) (see 


4.3 A note on the fitting strategy 

Our choice of separating the analysis in two steps (Secc. |4~T1 
and |4.2| ) is motivated by the fact that the structural param¬ 
eters (A7*,77*) physically describe the mass distribution of 
stellar discs and hence, in principle, they are best measured 
on the basis of available data for X* alone, irrespective of the 
distribution of newly born stars. Also, once such a measure¬ 
ment has been achieved, the fact that the SFRD profiles can 
be reproduced without a further tuning of (AT*, 77*) provides 
a valuable test for the validity of our theory. 

However, we also investigated whether our results would 
change if our 4 parameters (AT*, 77*, ^m, vb) were allowed to 
vary simultaneously to reproduce both the stellar mass and 
the SFRD radial profiles. To this purpose, we ran, for each 
galaxy, a Monte Carlo Markov Chain based on the com¬ 


bined likelihood of both our datasets (Secc. 3.2 and 3.31. 
We then compared the resulting radial growth rates with 
the ones derived with our preferred strategy, finding an ex¬ 
cellent agreement, with a median absolute difference of just 
2 • 10~ 4 GyP 1 . We found some discrepancy in just 4 cases, 
2 of which within 2a (NGC 1097 and NGC 7793) and the 
other 2 within 3cr (NGC 3184 and NGC 3351). Note that all 
these objects belong to the group for which an unweighted 
fit was found to provide a better description of the overall 


disc structure (see Sec. 4.1l, while the effect of oscillations 


induced by spiral structure was not taken into account in the 
MCMC experiment. This probably explains even the moder¬ 
ate discrepancies for this small subset. Furthermore, it shows 
that our partially subjective choice of the weights, discussed 
in Sec. |4.1| has a very limited impact on our results. 


To the first group belong distance and inclination. Dis¬ 
tance uncertainties affect the physical values of the derived 
mass and scalelengths, while the inclination uncertainty 
mainly affects the determination of the mass. Inclination 
also affects the normalization of both the stellar mass and 
SFRD profiles, but it does it exactly in the same way; it is 
easily seen that this implies a vanishing net effect on the es¬ 
timates of vm and vb, which are also, even more obviously, 
completely independent on the adopted distance. The im¬ 
portant consequence of this is that our method allows us to 
measure the specific mass and radial growth rates of discs 
with greater accuracy than the mass and scalelength them¬ 
selves, a fact that we will further exploit in Sec. [5] 

The second group of systematics comprises the mass-to- 
light ratio, the calibration of the FUV-to-SFRD conversion 
and the return fraction 1Z. Apart from second order effects, 
like possible variations of the mass-to-light ratio and the 
IMF with radius or morphological type, the main uncer¬ 
tainty coming from these systematics is a common multi¬ 
plicative factor for both the growth parameters, vm and vb, 
for the whole sample, or, in other words, a possible global 
rescaling of all the derived timescales. As examples of global 
systematics, we consider in some more detail the effect of 
the IMF and of the return fraction. In our calibrations, we 
have implicitly adopted a Salpeter IMF. To switch, for in¬ 
stance, to the more popular Kroupaj ( [200 1) IMF, we should 
divide the M/L ratio and hence all stellar mass surface den¬ 
sities by a factor 1.6 (see footnote [2j Sec. 3.21, while mul¬ 
tiplying all star formation rate surface densities by a factor 
0.63 ( Kennicutt fc Evans|2012 l. The net result on the sSFR 
(and hence on the estimates of vm and vb) is less than 1 %. 
This is due to the fact that in both cases the impact of the 
IMF is essentially driven by a common change in normaliza¬ 
tion associated to the contribution of very low mass stars. 
The effect of the return fraction 1Z is stronger: for instance, 
changing our adopted 7 Z — 0.3 into 1Z = 0.48 (which is the 
largest of the values suggested by Leitner fe Kravtsov|2011 1 
would imply a reduction of all growth rates by a factor 1.35 
and an equal increase of all timescales. Unfortunately, the 
return fraction is a quite uncertain parameter, since it is 
significantly affected not only by the IMF, but also by the 
details of the final-to-initial mass relation, which is very dif¬ 
ficult to determine observationally. However, we stress that 
the dimensionless ratio between vb and vm is unaffected by 
any of the systematics we have discussed so far. The impor¬ 
tance of this fact will be highlighted in Sec. [5] 


5 RESULTS 

The results of our analysis are listed in Table [2] The quoted 
errors are just the formal fitting ones; in particular, they 
do not take into account systematic uncertainties, which, as 
discussed in Sec. |4.4[ might be important for the structural 
parameters M* and 77*, but have a limited impact on the 
growth parameters v m and vb- 


4.4 Notes on systematics 

We distinguish between two kinds of systematics, those af¬ 
fecting individual galaxies in a different way and those af¬ 
fecting the whole sample more or less homogenously. 


5.1 Inside-out growth 

While our analysis is able to reveal both positive and nega¬ 
tive radial growth rates, we find that 32 galaxies, out of 35, 
show vb > 0. Of the remaining 3 galaxies with a formally 
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Table 2. Best fit structural (M* and R t ) and growth (i/m and i/r) parameters for galaxies in our sample. Formal fitting errors are 
reported, not including contributions due to distance, inclination and calibrations of conversion fomulae. Compared with stellar mass 
and scalelength, the growth parameters v m and i/r are less affected by systematic effects (see Sec:. |4.4[ i. Out of 35 studied galaxies, 32 
have a positive radial growth rate z/r. 


Galaxy 

M* 

(10 9 M©) 

R * 

(kpc) 

1 

S-H 

1 ° 

(M 

1 

o 

(10- 2 Gyr- 1 ) 

NGC 0024 

3.01 + 0.13 

1.62 ± 

0.02 

5.85 ± 1.03 

2.51 ± 0.52 

NGC 0337 

27.7 + 2.3 

2.15 ± 

0.06 

11.6 ± 3.5 

4.95 ± 1.72 

NGC 0628 

43.7+ 1.4 

3.64 ± 

0.05 

8.22 ± 0.36 

2.87 + 0.21 

NGC 0925 

11.3 ± 0.8 

3.97 + 

0.11 

10.3 + 0.5 

0.799 ± 0.381 

NGC 1097 

68.2 ± 8.1 

6.32 ± 

0.23 

8.74+ 1.04 

-2.05 ± 0.75 

NGC 1512 

14.7 + 2.2 

2.22 ± 

0.09 

3.80 ± 0.81 

1.22 ± 0.45 

NGC 1566 

78.0 + 5.3 

3.30 ± 

0.07 

8.21 ± 0.83 

2.90 ± 0.44 

NGC 2403 

7.19 ± 0.22 

1.51 ± 

0.02 

9.91 ± 0.25 

2.93 ± 0.16 

NGC 2841 

92.8 ± 3.4 

3.69 ± 

0.05 

1.62 ± 0.08 

0.612 ± 0.045 

NGC 2976 

2.25 ± 0.21 

0.802 ± 

0.028 

5.88 ± 0.57 

2.04 + 0.33 

NGC 3031 

48.3 + 2.5 

2.54 ± 

0.03 

1.99 ± 0.20 

0.750 ± 0.118 

NGC 3184 

17.2 ± 1.6 

2.42 ± 

0.08 

6.00 ± 0.80 

1.50 ± 0.55 

NGC 3198 

32.0+ 1.8 

3.65 ± 

0.07 

8.05 ± 0.81 

3.30 ± 0.42 

IC 2574 

1.21 ± 0.07 

3.01 ± 

0.08 

12.0 ± 1.1 

3.71 ± 0.70 

NGC 3351 

45.3 + 3.4 

2.86 ± 

0.05 

3.62 ± 0.17 

0.384 + 0.109 

NGC 3521 

54.2 ± 2.3 

1.85 ± 

0.03 

4.75 ± 0.56 

1.95 ± 0.29 

NGC 3621 

24.8 ± 1.0 

1.74 ± 

0.02 

9.84 + 2.05 

4.18 ± 1.04 

NGC 3627 

58.2 ± 2.5 

2.36 ± 

0.03 

4.66 ± 0.36 

1.24 + 0.21 

NGC 4236 

1.83 ± 0.10 

2.77 + 

0.07 

11.4 + 0.5 

4.70 ± 0.27 

NGC 4536 

27.0+ 1.7 

3.90 ± 

0.09 

11.2 + 0.4 

1.17 + 0.24 

NGC 4559 

42.7+ 1.4 

4.47 + 

0.05 

10.6+ 1.0 

4.10 ± 0.51 

NGC 4569 

78.3 + 3.3 

4.38 ± 

0.05 

2.07 + 0.11 

-0.02940 ± 0.0781 

NGC 4579 

87.6 ± 2.2 

3.56 ± 

0.03 

1.70 ± 0.24 

0.357 + 0.155 

NGC 4725 

117+ 11 

5.45 ± 

0.18 

2.17 + 0.25 

0.831 ± 0.155 

NGC 4736 

22.8 ± 3.1 

1.13 ± 

0.05 

4.75 ± 1.06 

1.18 ± 0.52 

NGC 4826 

51.0+ 1.3 

1.96 ± 

0.01 

1.55 ± 0.21 

-0.0850 + 0.1068 

NGC 5033 

23.5 ± 2.6 

3.91 ± 

0.15 

10.2 ± 0.7 

2.09 ± 0.42 

NGC 5055 

57.4 + 3.8 

2.50 ± 

0.05 

4.87 ± 0.63 

1.45 ± 0.34 

NGC 5194 

77.5 ± 6.1 

2.75 ± 

0.07 

7.56 ± 0.55 

1.26 ± 0.38 

NGC 5398 

4.86 ± 0.22 

1.95 ± 

0.03 

8.77 + 0.19 

3.57 + 0.10 

NGC 5713 

54.8 + 6.6 

1.94 ± 

0.06 

11.0 ± 3.6 

4.54+ 1.80 

IC 4710 

2.71 ± 0.22 

2.13 ± 

0.08 

7.61 ± 1.69 

1.97+ 1.04 

NGC 6946 

46.6 ± 2.2 

2.67 + 

0.05 

8.31 ± 0.54 

0.765 ± 0.529 

NGC 7331 

123 ± 13 

2.64 ± 

0.08 

10.5 + 3.8 

5.02 ± 1.88 

NGC 7793 

5.10 ± 0.23 

1.26 ± 

0.02 

9.10 ± 0.73 

2.58 ± 0.46 


negative radial growth rate, two (NGC 4569 and NGC 4826) 
have more than 100% uncertainty in i/r and hence are con¬ 
sistent with evolution of the stellar scalelength in one sense 
or the other, or with no evolution. Incidentally, we point 
out that both these galaxies are known to have peculiar 
properties: NGC 4569 is an anemic spiral in the Virgo clus¬ 
ter, probably significantly affected by ram pressure stripping 
( Boselli et al.||2006 i, while NGC 4826 is likely to have un¬ 
dergone a strongly misaligned merger, as suggested by the 
presence of a counter-rotating gaseous disc in the outskirts 


( Braun et al. 19941. For only one galaxy in our sample, NGC 


1097, we clearly find the signature of a shrinking of the disc. 
It is interesting to notice, a posteriori , that this galaxy has a 
very disturbed morphology. This is likely due to a strong in¬ 
teraction with the companion NGC 1097A[^] This object is 
listed as a peculiar elliptical in the RC3 catalogue (de Vau- 
couleurs et al.fl991), though with ‘uncertain’ classification. 


In the GALEX Atlas (Gil de Paz et al .|200 7|, NGC 1097A is 
clearly visible as a clump northwest of the prominent bar of 
NGC 1097, bright in NIR and NUV, faint in FUV and sur¬ 
rounded by an extended, FUV bright, disc-like structure, all 
properties common, in the GALEX Atlas, to the bulges of 
spiral galaxies. This is suggestive that the whole system may 
be a galaxy pair in an advanced stage of merging. We also 
noted that NGC 1097 is the object with the largest derived 
disc scalelength (6.32 kpc) and we verified that this result 
is not changed if we exclude from the exponential fit the 
whole radial range occupied by NGC 1097A. NGC 1097 has 
also been suggested to have undergone other significant in¬ 
teractions in the recent past (Higdon ik Wallin 20031. If our 
interpretation of a strong interaction state, likely a merger, 
for this system, is correct, then there is no surprise that 
it behaves differently from the regular evolution of isolated 
galaxies. 


3 Note that NGC1097A cannot be seen in our online Atlas, since 
it has been masked out from our 3.6 /im map. 


From these considerations, we can conclude that our 
findings are in excellent agreement with the general predic- 
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Figure 3. The specific mass and radial growth rates vm (top) and v\y (bottom) as a function of disc stellar mass (left) and morphological 
type (middle) and the relative histograms (right). Lower panels contain only the 32/35 galaxies with Vb > 0. Error bars are formal fitting 
uncertainties. The distributions of v m and vyt have some similarities (see text), but v\\ values are sistematically lower by ~ 0.5 dex. 


tions of the inside-out growth scenario for the evolution of 
spiral galaxies. 


5.2 Mass and radial growth rates 

In Fig. [3] the mass and radial growth rates of the galaxies 
in our sample are plotted against disc stellar mass and mor¬ 
phological type. Since we are using logarithmic units, radial 
growth rates (lower panels) are reported only for those 32/35 


galaxies with vb > 0 (see Sec. 5.11. Error bars represent 


formal fitting errors only. In particular, errors on distance 
and inclination are not taken into account in these plots. 
As discussed in Sec. |1.4[ these additional errors can affect 
stellar masses, but not vm and vb, which are only subject 
to a common multiplicative uncertainty due to global cali¬ 
bration issues. We also give in Fig. [3] the histograms for the 
distributions of vm and vb, binned in logarithmic intervals 
of 0.25 dex width. 

From the upper panels of Fig. [3] we see that the spe¬ 
cific mass growth rates (or specific star formation rates, or 
sSFR) of the discs of our galaxies take a relatively narrow 
range of values, with most of our points clustered around 
vm ~ 0.1 Gyr -1 , which corresponds to a mass growth 
timescale of ~ 10 Gyr. This is in substantial agreement 
with the typical sSFR of star-forming galaxies in the Local 
Universe (e.g. El baz et al.p Oll) and in particular with the 
relative constancy of the sSFR of the discs of spiral galaxies 
( Abramson et ah||2014 |, although a more detailed compar¬ 
ison would require more statistics and a careful treatment 


of global systematics (Sec. 4.41, which is beyond the scope 
of this work (see e.g. Speagle et al.|2014 about subtle issues 


concerning the homogeneization of measurements of these 
kind). We find that a small group of 6 galaxies (NGC 2841, 
NGC 3031, NGC 4569, NGC 4579, NGC 4826, NGC 4725) 
have a particularly low sSFR, with log(!z M /Gyr _1 ) < —1.5. 
We have already recognized two of them (NGC 4569 and 
NGC 4826) as objects with peculiar properties and a close 
to vanishing radial growth rate (see Sec. |5.1| ), but we cannot 
tell whether these peculiarities have a direct physical con¬ 
nection with the low measured values of vm■ However, we 
can see from the upper-mid panel that the whole group of 6 
slowly-evolving galaxies are also among the galaxies of the 
earliest types in our sample. This may be interpreted as an 
indication of downsizing (e.g. |Cowie et a l. 1 996| : galaxies 
with high mass and early types are more likely to have com¬ 
pleted most of their evolution in ancient epochs and hence 
to be growing with only mild rates nowadays. Also, galaxies 
of high mass and early-type might be more subject to star- 
formation quenching, the origin of which and its connection 
with morphology is still matter of investigation (e.g |Martig| 
et al.|2009 Pan et al.|2014 1 . 


In the lower panels of Fig. [3] we can see the distribution 
of the radial growth rates, which are the main novelty of this 
work. When plotted against disc mass, the radial growth rate 
vb shows a quite similar distribution with respect to the one 
of vm, but systematically shifted downwards by ~ 0.5 dex. 
This suggests that our galaxies are growing in size, on aver¬ 
age, at about 1/3 of the rate at which they are growing in 
stellar mass. The histogram of vb reveals a distribution that 
is similarly asymmetric, though less strongly peaked, with 
respect to the one of vm ■ More than 50% of our galaxies are 
in the two bins around log(;zR/Gyr -1 ) = —1.5, that is 0.5 
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Table 3. Basic statistics for our derived specific mass and radial 
growth rates (cfr. Fig. [3j. Note that discs with higher masses have 
lower median values and a higher scatter for both I'm and i'r. 



M* < 10 10 M© 

M* > 10 10 M© 

All 

log^M/Gyr x ) 

Median 

-1.05 

-1.12 

-1.09 

Scatter 

0.13 

0.25 

0.20 

l°g(i'R,/Gyr x ) 

Median 

-1.56 

-1.87 

-1.70 

Scatter 

0.18 

0.37 

0.35 


dex below the peak of the I'm distribution, corresponding 
to a typical radial growth timescale of ~ 30 Gyr. Since the 
radial growth rate of galaxy discs has been studied much 
less than the sSFR, it is less obvious to compare our find¬ 
ings with the ones of previous studies. However, we notice 
the typical timescale reported above is compatible with a 
radial growth of ~ 25% in the last ~ 7 Gyr, very similar to 
what found by Munoz-Mateos et al. (2011). Our typical i'r 
is also consistent with the results from the detailed study 


of resolved colour-magnitude diagrams for M 33 (Williams 


et al. 2009), while NGC 300 seems to have been growing 
at about half of this rate ( Gogarten et al.||2010 |. However, 
our estimates strictly refer to the current time and caution 
is mandatory in extrapolating these instantaneous measure¬ 
ments to a significant fraction of the past history of galaxies. 

Apart from the global vertical shift, the distributions 
of galaxies in the upper-left and bottom-left panels of Fig. 
[3] have a quite well defined common shape. In both cases, 
there is a horizontal upper envelope, close to the peak of the 
distribution, and a continuous increase of the scatter with 
increasing disc mass. Such a scatter is asymmetric and bi¬ 
ased towards low values of I'm and i'r, with the result that, 
on average, more massive discs appear to grow at a slower 
rate, both in mass and size, than the less massive ones. This 
effect is quantified in Table[3] where the median mass and ra¬ 
dial growth rates are reported, together with the associated 
scatter, for two subsamples with a derived disc mass lower 
or greater than 10 1(1 M©. It is not easy to understand, with 
our relatively small sample, if some simple physical property 
can be invoked to explain the scatter at high disc masses. 
However, we performed some simple checks and did not find 
any particular correlation between the position of galaxies 
in our plots and special properties, including the presence 
of a bar, an XUV disc, a break in the exponential profile, 
or indications of a warp or an interaction. We are therefore 
tempted to interpret the effect as intrinsic. 


6 IMPLICATIONS FOR THE EVOLUTION OF 
SCALING RELATIONS OF DISC GALAXIES 

6.1 The mass-radial growth connection 

Among our derived quantities, the mass and radial growth 
rates are the less affected by systematic uncertainties (see 


Sec. 4.41. Hence, the most reliable of our results are those 


that we can derive plotting I'm and i'r against each other. 
This is also important to understand whether the shift, that 
we found in Sec. 15.21 of a factor ~ 3 between I'm and i'r 



Figure 4. The relation between the specific mass growth rate 
and the specific radial growth rate of galaxy discs. The points are 
the results of our measurements, the lines are predictions of some 
simple models. The solid line is the expectation if known scaling 
relations of disc galaxies are not evolving with time, the dahed 
line and the dot-dashed line are for scaling relations evolving on 
a timescale of 100 Gyr, in one sense or the other (see text for 
details). A more rapid evolution is excluded by our results. 


is significant only at a statistical level or if it reflects an 
evolutionary property of individual galaxies. 

Indeed, this experiment (Fig. [4| reveals that the two 
growth rates are related to each other much more strongly 
than they are, individually, with mass or morphological type 
(cfr. Fig. [3}. 

The fact that masses and sizes of galaxies grow in an 
interlinked way is not very surprising on its own. Hence, 
rather than fitting a straight line to the points in Fig. [4j we 
prefer to seek for some simple physical explanation that can 
give a quantitative account to our finding. 


6.2 A comparison with a simple theoretical 
prediction 


Let us assume that mass and size of the discs of spiral galax¬ 
ies are connected by a power-law (e.g. |Courteau et al.|2007[ 
Lange et al.|2015 |: 

R* = AM: (8) 


Furthermore, let us assume that the coefficients A and a 
are not evolving with time, so that the relation Q defines 
not only the present locus, but also the evolutionary track 
of stellar discs. Then it immediately follows that the specific 
mass and radial growth rates should be linked by the very 
simple linear relation: 

i'r = a I'm (9) 

or, in logarithmic units: 


log i'r = log a + log I'm 


( 10 ) 


Independently on the value of a, equation (101 implies that 
(I'm , i'r) points should lie, in a double logarithmic plot like 
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Fig.g on a line of unitary slope; this is indeed the slope of 
the solid line drawn in Fig. [4] which gives a quite good ac¬ 
count for the distribution of our datapoints. This is already 
suggestive that our results are consistent with the existence 
of a non-evolving, power-law, mass-size relation for the discs 
of spiral galaxies. 

Of course, for our simple scenario to be fully predictive, 
not only the slope, but also the intercept, of such a straight 
line should be predicted as well, which is accomplished by 
specifying the expected value for a. To do this, we just com¬ 
bine two well-known scaling relations for disc galaxies, the 
Tully-Fisher relation (Mc. Gaugh|20121 , between the rotation 
velocity V and the mass M of a spiral galaxy: 


Uoc M u 


( 11 ) 


and the Fall relation (Romanowsky & Fall 2012), between 


specific angular momentum l and mass: 

l oc M o e 


( 12 ) 


Since exponential discs belong to a structurally self-similar 
family, one also has: 




(13) 


Substituting (121 and (111 into (131 we get a power-law 


mass-size relation of the form (|8|), with a = 0.35. This is 
not far from the value 0.32 empirically derived by |Courteau| 
et al. (2007) as an average slope for the mass-scalelength 


relation of disc galaxies in the Local Universe (since it was 
derived in the I-band, residual effects cannot be excluded 
arising from M/L variations). Shallower slopes are frequently 
found by studies based on half-light radius rather than disc 
scalelength (see e.g. |Lange et~a l.|[2015|; this is in qualita¬ 
tive agreement with expectations if we consider an obvious 
morphological effect (more massive galaxies tend to have 
more prominent bulges and hence a smaller ratio between 
half-light radius and scalelength). 

To define our simple model, we just retained the value 
a = 0.35, derived from the Tully-Fisher and the Fall rela¬ 
tions as explained above, and we adopted it to draw the solid 
line in Fig. [f] Hence we see that the majority of our data¬ 
points lie on a locus that can be independently predicted, 
without any free parameter, just assuming that known scal¬ 
ing relations for disc galaxies hold and are not evolving with 
time. Since these simple hypotheses are completely indepen¬ 
dent from the way our results were derived, the agreement 
between the two is very unlikely to occur by chance, or to 
be due to biases of any kind, and we are tempted to inter¬ 
pret this finding as an indication for the validity of both our 
method and the hypotheses themselves. 

We can also consider what effect residual systematics on 
vm and un could have on our findings. As discussed in Sec. 
|4.4| the effect is an unknown common multiplicative factor 
for i/m and v^. In the diagram shown in Fig. [4] this implies 
a collective motion of all points along a line of unitary slope, 
or, equivalently, a mapping of the theoretical line into itself. 
Therefore, our conclusions are robust at least against the 
most obvious systematic uncertainties. 

A word of caution is appropriate, however, against pos¬ 
sible overinterpretation of our model and result. In fact, the 
Tully-Fisher relation is known to hold better for the whole 
baryonic content of spiral galaxies (McGaughp~012 I, while 


we have impicitly applied it just to the stellar mass of the 
disc. On the other side, the Fall relation seems to hold better 
when the disc component is considered separately from the 
bulge. Hence, it can be argued that, in deriving our predic¬ 
tions, we have mixed non-homogeneous empirical evidence. 
A more detailed analysis, taking this aspect into a proper 
account, would be interesting, but is beyond the scope of 
this work. 


6.3 Evolutionary effects 


In the previous Section we have seen that our results are 
compatible with the Tully-Fisher and Fall relations to be 
not evolving with time. To quantify this statement, we put 
here an upper limit on how fast a possible evolution can be in 
order to be still compatible with our results. For simplicity, 
we focus our attention on the evolution in normalization, 
although a similar analysis could be performed for the slope 
evolution as well. 

If, in ( |8|) , we allow the normalization A to change with 
time, then (|9|) simply modifies into: 


where: 


= va + atv M 


v A (t) := |(ln A)(t) = 


(14) 


(15) 


is the specific growth rate of the normalization coefficient 


A. Of course, when equation (141 is compared with our ob¬ 


servations in the Local Universe, va has to be intended as 
evaluated at the present time. 

The dashed and dot-dashed lines in Fig. [4] show the 
predictions of two models with a very mild evolution in nor¬ 
malization, in one sense or the other: va — + 0.01 Gyr -1 . It 
is clearly seen that the predicted distribution of galaxies in 
the (i'm,J / r) plane is extremely sensitive to the parameter 
va , making this diagram a new and powerful observational 
tool to constrain the evolution of scaling relations of galaxy 
discs. Also, since the two additional lines are both inconsis¬ 
tent with the empirical distribution, we quantitatively infer 
that, even admitting that scaling relations are evolving with 
time, they are doing so on timescales that are larger than 
(v A , max) — 1 = 100 Gyr, hence much larger than the Hubble 
time. 

Strictly speaking, the statement above mainly refers to 
the mass-size relation (|8|. In fact, although the Tully-Fisher 


relation (111 and the Fall relation (121 are the backbone of 


the simple model sketched in |6.2| it may be considered not 
trivial to draw out conclusions concerning them individually, 
since they both involve kinematics, while we did not directly 
make use of kinematical data. However, our results indicate 
that an evolution of the Tully-Fisher relation, if present, 
has to be accompanied and finely balanced by an opposite 
evolution of the Fall relation. 

If compared with the direct observational study of scal¬ 
ing relations of disc galaxies at different redshifts (with 
all the appropriate caveats, see Sec. [I]), our results are in 
better agreement with those finding little or no evolution 
(see again references in Sec. [I]), for either the mass-size or 
the Tully-Fisher relation, while, to our knowledge, no sim¬ 
ilar studies are available yet concerning the Fall relation. 
However, we stress again that our empirical upper limit 
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[\va\ < 0.01 Gyr -1 ) only refers to v a evaluated at the 
present time and hence it is related (though not equal) to 
the slope of empirical A(z) relations as measured at z = 0. 
For instance, our results have no formal tension with those 
of Trujillo et al. (2006), who still find a preference for evo¬ 
lutionary models, since their data points are also perfectly 
consistent with the slope of the A(z) relation to vanish up 
to z ~ 1; to tell the difference between models, more precise 
measurements at moderate redshift, or the use of a sensitive 
local diagnostics like the one that we have proposed here, 
can be a valuable complement to pioneering observational 
campaigns in the extremely distant Universe. 


7 SUMMARY 

In this work, we have developed, from very simple assump¬ 
tions, a model that predicts a universal shape for the radial 
profile of the star formation rate surface density (SFRD) of 
spiral galaxies. This model accounts for the basic properties 
of observed profiles and naturally includes a parametrization 
of the growth of stellar discs. As a consequence, we have de¬ 
vised a novel, simple and powerful method to measure the 
instantaneous mass and radial growth rates of stellar discs, 
based on their SFRD profiles. We have applied our method 
to a sample of 35 nearby spiral galaxies. Our main results 
are: 

(i) For most of the galaxies in our sample, the SFRD pro¬ 
file is satisfactorily reproduced by our model, in such a way 
that we could measure the mass and radial growth rates I'm 
and I'd, of their stellar discs. 

(ii) Virtually all galaxies show the signature of inside-out 
growth (i/r > 0 ). 

(iii) Typical timescales for the mass and radial growth of 
our stellar discs are of the order of ~ 10 Gyr and ~ 30 Gyr, 
respectively, with some uncertainty due to systematic ef¬ 
fects. 

(iv) The mass and radial growth rates appear to obey a 
simple linear relation, with galaxy discs growing in size at 
~ 0.35 times the rate at which they grow in mass. Compared 
with the individual timescales given above, this dimension¬ 
less ratio is more robust against systematic uncertainties. 

(v) The distribution of galaxies in the (I'm, i'r) plane is a 
sensitive diagnostics for the evolution of scaling relations of 
galaxy discs. 

(vi) Our results are in very good agreement with a simple 
model, without free parameters, based on the universality 
of the Tully-Fisher relation and the Fall relation, suggesting 
that they are not evolving with time. Possible residual evo¬ 
lution is constrained to occur on timescales that are much 
larger than the age of the Universe. 
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